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ABSTRACT 

We propose a new version of the Baade-Becker- Wesselink technique, which allows 
one to independently determine the colour excess and the intrinsic colour of a radially 
pulsating star, in addition to its radius, luminosity, and distance. It is considered 
to be a generalization of the Balona approach. The method also allows the function 
F(CIq) = BC(CIq) + 10 x log(T e ff(CIo)) for the class of pulsating stars considered 
to be calibrated. We apply this technique to a number of classical Cepheids with 
very accurate light and radial-velocity curves and with bona fide membership in open 
clusters (SZ Tau, CF Cas, U Sgr, DL Cas, GY Sge), and find the results to agree well 
with the reddening estimates of the host open clusters. The new technique can also 
be applied to other pulsating variables, e.g. RR Lyraes. 
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1 INTRODUCTION 

Classical Cepheids are the key standard candles, which are 
used to set the zero po int of the extragalactic distance scale 
jFreedman et al.lboOll ') and a lso serve as young-populatio n 
tracers of great importance (Bi nnev and Merrifieldl fl998). 
They owe their popularity to their high luminosities and 
photometric variability (which make them easy to identify 
and observe even at large distances) and the fact that the 
luminosities, intrinsic colours, and ages of these stars are 
closely related to such an easy to determine quantity as the 
variability period. 

It would be best to calibrate the Cepheid period- 
luminosity (PL), period-colour (PC), and period- luminosity- 
colour relations via distances based on trigonometric par- 
allaxes, however, the most precisely measured parallaxes 
of even the nearest Cepeheids remain insufficiently accu- 
rate and, more importantly, they may be fraught with so 
far uncovered syste matic errors. Here the Baade-Becker- 
Wess elink method (jBaadd 1 19261 ; iBeckeij Il940l : IWesselinkl 
1 19461 ) comes in handy, because it allows the Cepheid dis- 
tances (along with the physical parameters of these stars) 
to be inferred, thereby providing an independent check for 
the results based on geometric methods (e.g., trigonometric 
and statistical parallax). 

However, all the so far proposed versions of the 
Baade -Becker- Wesselink meth od (surface brightness tech- 
nique jBarnes and Evanslll976h . maximum-likelihood tech- 



nique jBalondirtmr n depend, in one way or another, on 
the adopted reddening value. Both techniques are based on 
the same astrophysical background but make use somewhat 
different calibrations (limb-darkened surface brightness pa- 
rameter, bolometric correction - effective temperature pair) 
on the normal co lour. Here we propose a generalization of 
the iBalonal (|l9T7t ) technique, which allows one to indepen- 
dently determine not only the star's distance and physical 
parameters, but also the amount of interstellar reddening, 
and even calibrate the dependence of a linear combination 
of the bolometric correction and effective temperature on 
intrinsic colour. 



2 THEORETICAL BACKGROUND 

We now briefly outline the method. First, the bolometric 
luminosity of a star at any time instant is given by the fol- 
lowing relation, which immediately follows from the Stefan- 
Boltzmann law: 



L/L Q = (R/Rq) 2 x (T/T @ ) 4 



(1) 



Here L, R, and T are the star's bolometric luminosity, ra- 
dius, and effective temperature, respectively, and the sub- 
script denotes the corresponding solar values. Given that the 
bolometric absolute magnitude Mbol is related to bolometric 
luminosity as 



M bo i = M bo iQ - 2.5 x log(L/L e ), 



(2) 
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we can simply derive from Eq. (1): 
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M bol - M bolQ = -5 X log(R/R e ) - 10 x log(T/T Q ) (3) 

Now, the bolometric absolute magnitude Mbol can be writ- 
ten in terms of the absolute magnitude M in some pho- 
tometric band and the corresponding bolometric correction 
BC: 

M bol = M + BC, (4) 

and the absolute magnitude M can be written as: 

M = m - A - 5 x log(d/10 pc). (5) 

Here m, A, and d are the star's apparent magnitude and in- 
terstellar extinction in the corresponding photometric band, 
respectively, and d is the heliocentric distance of the star in 
pc. We can therefore rewrite Eq. (3) as follows: 

m = A + 5 x log (d/ 10 pc) + Mbol® + 10 x log(T Q ) 

- 5 x log(R/R e ) - BC - 10 x log(T). (6) 

Let us introduce the function F(CI ) = BC + 10 x log(T), 
the apparent distance modulus (m — M) app = A + 5 x 
log(d/10 pc), and rewrite Eq. (7) as the light curve model: 

m = Y - 5 x log{R/R Q ) - F. (7) 

where constant 

Y = (m - M)a PP + Mboie + 10 x log(T e ). 

We now recall that interstellar extinction A can be de- 
termined from the colour excess CE as A = R\ x CE, 
where R\ is the total-to-selective extinction ratio for the 
passband-colour pair considered, whereas M bo io, Rq, and 
Tq are rather precisely known quantities. The quantity 
F(CIo) = BC + 10 x log CT) is a f unctio n of intrinsic colour 
index CI = CI - CE. iBalonal (|l977h used a very crude 
approximation for the effective temperature and bolomet- 
ric correction, reducing the right-hand of the light curve 
model (Eq. 7) to the linear function of the observed colour, 
with the coefficients containi ng the colour excess i n a latent 
form. It should be noted that lSachkov et ail (1998); Sachkov 
(2002) used non-linear approximation in Eq. (7) to calculate 
Cepheid radii. 

The key point of our approach is that the values of func- 
tion F are computed from the already available calibrations 
of bolometric correction BC(CIq) and effec t ive tempera- 
ture T(CIn) (lFlowerlll996l:lBessell et al.lll998l;lAlonso et al 



1999; Sckiguchi and Fukugita 2000; Ramirez and Melcndcz 
2005: [Biazzo et al.ll2007l ;l Gonzalez Hernandez and Bonifacid 



20091 ). These calibrations are expressed as high-order power 
series in the intrinsic colour: 



F{CI ) = a + ^a k -CI%, 



(8) 



with known {ofe} and N amounting to 7; in some cases the 
decomposition also includes the metallicity ([Fe/H]) and/or 
gravity (log g) terms. 

As for the star's radius R, its current value can be de- 
termined by integrating the star's radial- velocity curve over 
time (dt = (P/2tt) ■ dip): 



where Ro is the radius value at the phase <po (we use mean 
radius, < R >= (R m in + Rmax)/2 ); Vy, the systemic radial 
velocity; <p, the current phase of the radial velocity curve; P, 
the star's pulsation period, and pf is the projection factor 
that accounts for the difference between the pulsation and 
radial velocities. Given the observables (light curve - appar- 
ent magnitudes m, colour curve - apparent colour indices 
CI, and radial velocity curve - V r ) and known quantities 
for the Sun, we end up with the following unknowns: dis- 
tance d, mean radius < R >, and colour excess CE, which 
can be simply found by the least-squares or maximum like- 
lihood technique. 

In the case of Cepheids with large amplitudes of light 
and colour curves (AC/ ^ 0.4 m ) it is also possible to apply a 
more general technique by setting the expansion coefficients 
{at} in Eq. (8) free and treating them as unknowns. We 
expanded the function F = BC + 10 x log(T) in Eq. (7) 
into a power series about the intrinsic colour index CIq 1 of 
a well-studied "standard" star (e.g., a Per or some other 
bright star) with accurately known T st : 



F = BC 3t + 10 x log(T s 



N 

+ ^2a k -(CI 



CE- 



■C7o') fe (10) 



The best fit to the light curve is provided with the 
optimal expansion order TV ~ 5 — 9. We use this modifi- 
cation to calculate the physical parameters and reddening 
CE of the Cepheid, as well as the calibration F(CIo) = 
BC(CIo) + 10 x log(T e ff(CIo)) for the given metallicity 
[Fe/H]. 



3 OBSERVATIONAL DATA, CONSTANTS, 
AND CALIBRATIONS 

Our sources of data include Berdnikov's extensive mul- 
ticolor photoel ectric and CCD photo metry of classi- 
cal Cepheids l|Berdnikovl 1 19951 . 120081 ) and very ac- 
curate ra dial-velocity measu r emen t s of 165 northern 
Cepheids (|Gorvnva et al.lll992l . [l99^ . Il998l . |2002J ) taken in 
1987-2009 (about 10500 indi vidual observati ons) with a 
CORAVEL-type spectrometer (jTokovininlll987f ). These data 
sets are nearly synchronous, to prevent any systematic errors 
in the computed radii and other parameters due to the evo- 
lutionary period changes resulting in phase shifts between 
light , colour and radial velocity variations (I Sachkov et al. 
199Sh . We adopt T Q = 5777 K, M bole = +4.76 m (|GraT 



R(t) -Ro = -pf ■ 



(V r (t) - V y ) ■ (P/27T) • dip, 



(9) 



2005J). We proceeded from (V, B — V) data and found the 
best solutions for the V-band light curve and (73 — V) color 
curve to be those computed usi ng the F((B — V)p) func- 
tion based on two calibrations jFlowej 1 199^ :1 Bessell et al.l 
1998) of similar slope (see Fig. 2 e); the poorer results ob- 
tained using the other cited calibrations can be explained 
by the fact that the latter involved insufficient number of 
supergiant stars. 



4 THE PROJECTION FACTOR 

There is yet no consensus concerning the p rojection factor 
(PF) value to be used for Cepheid variable s (Nardctt o et al.l 
12004 lGroenewegenll2007l : iNardetto et al.ll2007l . 120091 ). Dif- 
ferent authors use constant values ranging from 1.27 to 
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1.5, as well as PFs depending on the pulsation period and 
other parameters. Different approaches lead to small sys- 
tematic differences in the inferred Cepheid parameters, first 
of all, in the rad ii. Based on geometrical considerations, 
iRastorguevI |2010) derived phase-dependent PFs as sim- 
ple three-parametric analytic expressions depending on the 
pulsation velocity, limb darkening coefficient, and spectral 
line broadening, adjusted to CORAVEL radial velocities 
of Cepheids. We suspect that the period dependence re- 
flects mainly the dependence of the PF on limb darkening. 
To compare our results with other calculations, we finally 
adopted a moderate dependen ce of PF on the period advo- 
cated bv lNardetto et all \200it i: 

p = (-0.064 ± 0.020) x log{P, days) + (1.376 ± 0.023), (11) 

though we repeated all calculations with other variants of 
PF dependence on the period and pulsation phase to assure 
the stability of the calculated colour excess. 



5 PRELIMINARY RESULTS 

To test the new method, we used the maximum likelihood 
technique to solve Eq. (7) for the V-band light curve and 
B — V colour curve for several classical Cepheids resid- 
ing in young open clusters: SZ Tau (NGC 1647), CF Cas 
(NGC 7790), U Sgr (IC 4725) PL Cas (NG C 129), GY Sge 
(anonymous OB-association (Forbes 1982)) as well as for 
approximately 30 field Cepheids from o ur samp le . We found 
two log(T e ft) cali brations - those of iFlowerl (|l996ft and 
iBessell et all (|l998h - combined with the BC(V) cali bration 
as a fu nction of normal colour (B — V)o proposed by IFlowerl 
( 1996) - to yield the best fit to the observed V-band light 
curve via Eq. (7). A weak sensitivity of calculated reddening, 
Eb-v, to the adopted PF value (constant or period/phase- 
dependent) and to the derived < R > value can be explained 
by very strong dependence of the light curve's amplitude on 
the effective temperature, ~ 10 x log(T), and, as a conse- 
quence, on the dereddened colour. 

Though the internal errors of the reddening Eb-v seem 
to be very s mall, the value s dete rmined using t he tw o best 
calibrations, IFlowerl (| 1996ft and IBessell et all l|l998ft . may 
differ by as much as 0.03 — 0.05 m , due to the systematic 
shift between these two calibrations (Fig. 2 e). Table [T]sum- 
marizes the inferred parameters for the cluster Cepheids 
studied. Fig. 1 shows the observed and smoothed data and 
the final fit to the V-band light curve for U Sgr Cepheid. 
Our reddenings seem to agree well with the corresponding 
WEBDA values, particularly if we remember that the er- 
rors of the adopted cluster reddening estimates are as high 
as ±0.05™. Our next step will be to make use of the cali- 
brations of T e f f and BC as a function of red and infrared 
colours (V — R, V — I, V — K) and to compare derived 
reddening ratios with the conventional extinction laws. 

Note that the inferred radius and luminosity of SZ Tau 
are too large for its short period; this Cepheid probably pul- 
sates in the 1 st or even in the 2 nd overtone, as may be indi- 
rectly evidenced by its low colour amplitude (about 0.15 m ). 

Fig. 2 shows the observed data, the fit to the V- 
band light curve, and the inferred calibration F = 10 x 
log(T eff ) + BC{V) vs (B - V)o) calculated for TT Aql 
Cepheid (as a 5*' 1 -order expansion in the normal colour). 



The i nferred calibration is very close to that of IFlowerl 
(|l996h . We used a Per as the "standard" star, with T 3t » 
(6240 ± 20) K, [Fe/H] « -0.28 ± 0.06 i|Lee et alj|2006ft . 
(B - V) st « 0.48 m and E B -v » 0.09 m (WEBDA, for 
a Per cluster). To take into account the effect of metal- 
licity on the zero-point F(CIo) st , we estimated the gra- 
die nt dF(CI ) st I d YFeJH ] w +0.24 from the cali b ration s 
by lAlonso et al.1 i 19991) ; ISekiguchi and Fukugital l|2000h ; 
iGonzalez Hernandez and Bonifacio! l|2009ft . For TT Aql, 
Eb-v ~ (0.65 ± 0.03) m . In some cases (with large ampli- 
tude of color variation) the "free" calibration (Eq. 10) can 
markedly improve the model fit to the observed light curve 
of the Cepheid variable. Fig. 2 f shows the example of cali- 
brations of the F functions derived from nine Cepheids with 
different metallicity and surface gravity values. Temperature 
difference at T e ff ~ 6600 — 5100 K is amounted to 3 — 5%. 

When applied to an extensive sample of Cepheid vari- 
ables with homogeneous photometric data and detailed ra- 
dial velocity curves, the new method is expected to give a 
completely independent scale of reddenings, a new Period - 
Colour - Luminosity relation, and a new distance scale for 
the Milky- Way Cepheids. 
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Table 1. Physical parameters, distances, and interstellar reddening values for the cluster Cepheids 
analysed using the new version of the BW method. Reddening values from WEBDA dat a base 
{http://www.univie.ac.at/webda/l l are also shown for comparison; asterisk: Eg_y adopted from I Forbes! 
{1983). Distances are calculated using Ry = 3.3. 



Star 


Cluster 


Period (d) 


Distance (pc) 


E B -v 


E B - V (WEBDA) 


< R > /R Q 


M v 


SZ Tau 


NGC 1647 


3.149 


796±90 


0.40±0.02 


0.370 


57.0±7.0 


-4.32±0.25 


CF Cas 


NGC 7790 


4.875 


3585±87 


0.54±0.02 


0.531 


46.7±0.9 


-3.41±0.05 


U Sgr 


IC 4725 


6.745 


613±25 


0.50±0.03 


0.475 


54.2±1.8 


-3.90±0.08 


DL Cas 


NGC 129 


8.001 


2067±58 


0.47±0.05 


0.548 


69.3±1.6 


-4.12±0.06 


GY Sge 


Anon OB 


51.78 


2136±163 


1.44±0.05 


1.29±0.06 (*) 


208±11 


-6.27±0.15 




Figure 1. Panel (a): Observed and fitted radial-velocity curve of U Sgr. Standard deviation uy r = 1.1 km/s. Panel (b): Observed and 
smoothed colour curve. Panel (c): Radius variation with phase. Panel (d): Observed and fitted light curve. 
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Figure 2. Panel (a): Observed and fitted radial-velocity curve of TT Aql. Standard deviation ay r = 1.3 km/s. Panel (b): Observed 
and smoothed colour curve. Panel (c): Radius variation with phase. Panel (d): Observed and fitted light cu rve. Panel (e): Calculate d 
calibration for TT Aql (function F = 10 X log(T eff ) + BC(V) vs (B - V) )) and calibrations by |Flowej i | 19961) and lBessell et al.l GUI). 
Also shown is the position of the standard star a Per corrected for metallicity difference. Panel (f): Calculated calibration (function 
F = 10 X log(T e ff) + BC(V) vs (B — V)o)) for 9 Cepheids with large amplitudes of the colour curves and different metallicities. 



